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We present a theoretical estimation for the cross-section of exclusive two p-meson production in two 
photon collision when one of the initial photons is highly virtual. The compatibility of our analysis 
with recent experimental data obtained by the L3 Collaboration at LEP is discussed. We show 
that these data prove the scaling behaviour of the exclusive production amplitude. They are thus 
consistent with a partonic description of the exclusive process 77* — > p°p° when Q 2 ~ 2 — 20 GeV 2 . 



INTRODUCTION 

Two-photon collisions provide a tool to study a variety of fundamental aspects of QCD and have long been 
a subject of great interest (cf., e.g., 

QBB 

and references therein). A peculiar facet of this interesting 
domain is exclusive two hadron production in the region where one initial photon is highly virtual (its virtuality 
being denoted as Q 2 ) but the overall energy (or invariant mass of the two hadrons) is small |^|. This process 
factorizes 

BB 

into a perturbatively calculable, short-distance dominated scattering 7*7 — > qq or 7*7 — > gg, 
and non-perturbative matrix elements measuring the transitions qq — + AB and gg — * AB. These matrix 
elements have been called generalized distribution amplitudes (GDAs) to emphasize their close connection to 
the distribution amplitudes introduced many years ago in the QCD description of exclusive hard processes Q- 
In this paper, we focus on the process 7*7 — > p°p° which has not been much discussed theoretically (see however 
Ref. 0) but has recently been observed at LEP in the right kinematical domain Q. 

Since 7*7 — > pp is the crossed channel of virtual Compton scattering on a spin-1 meson, the physics considered 
here is closely related to deeply virtual Compton scattering (DVCS) on a spin-1 target [l{J, which has recently 
attracted some attention in the context of generalized (skewed) parton distributions of the deuteron . 
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KINEMATICS 
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The reaction which we study here is (see, Fig. 1): 

e(k) + e(Z) -> e(fc') + e(0 + P °( Pl ) + p°(p 2 ) (1) 

where the initial electron e(k) radiates a hard virtual photon with momentum q — k — k' , in other words, the 
square of virtual photon momentum q 2 = —Q 2 is very large. This means that the scattered electron e(fc') is 
tagged. To describe reaction (JTJ, it is useful to consider, at the same time, the sub-process : 

e{k)+ 1 { q ')^e{k') + p°{p l )+p (p 2 ). (2) 

Regarding the other photon momentum q' = I — I', we assume that, firstly, its momentum is collinear to the 
electron momentum I and, secondly, that q' 2 is approximately equal to zero, which is a usual approximation 
when the second lepton is untagged. 

Let us now pass to a short discussion of kinematics in the 77* center of mass system. We adopt the z axis 
directed along the three-dimensional vector q, and the p-meson momenta lie in the (x, z)-plane. It means that 
we ignore the azimuthal dependence of the cross-section, which happens to be absent in the approximation we 
use. So, we write for the momenta in the cm. system : 

q = (g , 0, 0, q), pi = (pi, pisintf, 0, picosfl). (3) 

Also, we need to write down the Mandelstam ^-variables for the electron-positron (|TJ and electron-photon (J2J 
collisions: 

S ee = (k + l) 2 , S el = {k + q'f. (4) 
Neglecting the lepton masses, these variables can be rewritten as 

S ee w 2(k ■ I), S el w 2(k ■ q') = x 2 S ee , (5) 
where the fraction x 2 defined as q' Q — x 2 Iq is introduced (see, |12|). 




Figure 1: Kinematics of the process e(k) + e(Z) — > e(fc') + e(Z') + p°(pi) + p°{p2) in the c.m.s of the two mesons. 



PARAMETERIZATION OF p-MATRIX ELEMENTS AND THEIR PROPERTIES 
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Let us first introduce the basis light-cone vectors. We adopt a basis consisting of two light-like vectors p and n 
of mass dimension 1 and — 1, respectively. In other words, they obey the following conditions : 

p 2 =n 2 =0, (p-n) = l. (6) 

With the help of this basis, the p-mesons momenta pi and P2 can be written as 

W 2 A£ 
^ = ^ + (1-0^-^, 

W 2 A M 

^ = (l-0/ + Cy^ + f (7) 

As usually for the two meson generalized distribution amplitude case, we introduce the sum and difference of 
hadronic momenta which take the form in the light-cone decomposition : 

A"=^-pS' ) P"=^+pS*. (8) 
The skcwcdncss parameter ( is defined by 



p+ _ l + ffcosfl / Am 2 p (q] 

Note also that within the cm. frame the transverse component of the transfer momentum At = (0, At, 0) is 
given by 

A T = -A T = [ Am 2 - W 2 ) sin 2 (9. (10) 



Let us now write down the decomposition for the longitudinal and transverse p-meson polarization vectors: 

(0) _ 1 {„ m P„ 



(i) 2 /V ^ C(2C"1)^ 2 + A^ \ 
- 1 CP2 M -(1-C)Pia« ^7 n M h 



(2) 2 /n\ 

6 lM = 7=^e M "/37P2aPl/3n 7 , (11) 

V T 



for the p-meson with momentum pi , and 



,(0)_ 1 f„ m P 



(i) 2 ^ , (\-Q(X-\)W 2 - A 2 T 

4 1 = ( (i - Opim - CP2, + ^ — 2(i-c) " 



e 



(2) 2 



= e M a/3 7 P2aPl/3« 7 , (12) 



-2 /j 

V "T 

for the other p- meson with momentum p 2 . As usual, 

ei- Pi = 0, ■ e\ X,) = -6 xy . (13) 



Besides, for each meson, the following polarization vectors can be introduced to define the light-cone helicity f 

e(0) = e(°), e(±) = T£ ^ . (14) 

We now come to the parameterization of the relevant matrix elements. Keeping the terms of leading twist 2, 
the vector and axial correlators can be written as: 

(pi,Ai;j>2,A a |^(0)7^(An)|0) = e?elV^( Pl , P2 ,n)HeP' V (yX,W 2 ), (15) 

i 

(p 1 ,A 1 ;p 2 ,A 2 |^(0) 7M 75^(An)|0) £ P p £ e^e p 2 A%{ PllP2l n)H? p > A (y, £ W^ 2 ). (16) 

Here Ai and A 2 are the helicities of p mesons and = denotes the Fourier transformation with measure {z\ = 
Xn,z 2 = 0) 0: 

d^{y) = dye- %ypzl+i ^ v)pZ2 . (17) 

With the help of parity invariance we can show that the vector tensors V^l may be written in terms of five 
tensor structures while the axial tensors A^l are linear combinations of four independent structures. In complete 
analogy with the analysis of the deuteron generalized parton distributions |ll| we write 

E e ? e ^i?(Pi:f2,n)ifr v (y) = -(e, • e 2 ) \y) + 

i 

(ei • n)(e 2 • Pl ) + (e 2 • n)( ei ■ p 2 )) H? > V (y) - ^ ' ' Pl) H^ v (y) + 

((ei • »)(ea • Pl ) - (e 2 • n)( ei .p 2 )) fl™ y (j,) + Urn 2 ( ei • n)(e 2 • n) + i(e x • e 2 )) ^ V (y) (18) 
for the vector tensor structures and 

J2^2 A %(P^P2^)Hr A (y) = -iel^ a e^Hr A (y) + 

i 

■ t * eiflfo -pi) +e 2/3 (ei -p 2 ) A • t A e^(e 2 ■ Pl ) - e 2p {e 1 ■ P2 ) A 
H nip ' nip 



ie£pA a (eip(e 2 • n) + e 2 p{e 1 ■ n) \ H% p ' A (y) 



(19) 



for the axial tensor structures. In 118J1 and (|19|l . the standard notation e^g = e a p~ ( s P ^ns was used and the 
dependence of parameterizing functions (GDA) on £ and W 2 is implied. 

Let us now turn to the consideration of symmetry properties. The 77* subprocess is selecting the parts with 
the following symmetry: 

Hr v (y) = -Hr v (i-y), Hr A (y) = Hr A (i-y) (20) 

Note that the scale dependence of generalized distribution amplitudes acquired in the process of factorization of 
the scattering amplitude has been studied in (I2j and there are no essential differences between the ttit channel 
discussed there and vector contributions to our pp case. At the same time, the evolution of axial contributions 
are similar to the case of the distribution amplitudes of singlet axial mesons. The W 2 behaviour of 27r GDA's 
has been recently explored in terms of an impact representation of the hadronization process qq — > MM jl^ . 



QUARK-HADRON HELICITY AMPLITUDES 
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Although the main objects of our investigation are the generalized distribution amplitudes (GDA) let us begin 
from the discussion of the helicity amplitudes related to the generalized parton distribution (GPD), which are 
related to GDAs by means of s t crossing symmetry ^[ |2 • One considers the quark helicity conserving 
distributions parameterizing the combination of the vector and axial matrix elements : 



A 



(A 2 ± ;Ai±) 



-(p 2 ,A 2 |V;(Q)7 + (l±75Mz)bi,Ai 



(21) 



where A and fi — ± denote the helicities of initial (final) hadrons and quarks, respectively. The matrix elements 
in l|21() are taken between two p mesons. In the forward limit when pi = p 2 helicity conservation takes place 
and requires Ai + /ii = A2 + ^2- The parity transformation, Ai — > — Xi etc., invariance leads to nine independent 
helicity amplitudes 1)21(1 . At the same time the time reversal transformation, (pi, Ai) <-> (/12, A2), invariance does 
not lead to any reduction of the number of independent structures owing to £-dependence (see, for instance, 

Q>. 

As noted before, the crossing transformation relates the helicity amplitudes referring to the GDAs to the 
corresponding GPDs helicity amplitudes. Indeed, the crossing transformations imply that the initial hadron 
is replaced by a final hadron with the opposite momentum : pi — > —pi and, therefore, opposite helicities 
Ai — > — Ai. Similarly, we implement the crossing replacement for the quark fields. Namely, the final quark with 
helicity fi2 is replaced by the initial antiquark with helicity — /i2- Thus, the quark helicity non-flip amplitudes 
in i-channel come to the amplitudes in s-channel where the initial quark and anti-quark helicities are opposite, 
and vice versa : 



A {t) 



A 



(») 

(A 2 -Ai ;T±)' 



(22) 



Note that the first two indices labeling the helicity amplitudes in IL'l'l) or (|21|) correspond to the helicities in the 
final state and the last two indices - to the initial state. 

Let us now focus on the amplitudes in s-channel. As the helicity and chirality of antiquarks are distinguished by 
sign (in this paper, we consider the case of massless quarks), the chirality conserving quark-antiquark operator 
will give the combination of quark and antiquark fields with opposite helicities : u+7 + u+ <-> «(_)7 + it( + ), where 
the bracketed subscripts denote the helicity while the unbracketed ones denote the chirality. Therefore, using 
(|15|) and l(16|l . we write, forgetting from now on the (s) subscript, 



A(x 2Xl ; ± T) = \{P2, X 2 ;pi, A 1 |^(0) 7 + (1 ± 7 5 )V^)|0). 



(23) 



Further, a straightforward calculation derives the expressions for the helicity amplitudes (cf. , ) 



2A++.+-) = H P i P ' V 



BH. 



pp,v 



PP,V 



Hi 



pp,A 



2B 



H. 



p P ,A 



(2C - 1)H. 



pp,A 
3 



2A 



(0+,+-) 



2BC i H r v -c 



H. 



pp,v 



H 



pp,v 



(2C- 1)W 2 



Am 2 p 



+ 2(1-0 



H. 



pp,v 



-H, 



pp,v 



+ (1 - CX"' 4 + (1 - 



(2C - l)W 2 



+ 4(1-C)£ 



H: 



pp,A 



H. 



pp,A 



2A 



(-+,+-) 
2^4(oo.+-) 



C(i - CM 



pp,A 



2B[ -H ppV -(2{-l)K 



PP' A — jjPP> A 



1 \ W 2 



H[ 



pp,v 



1 

-H 



pp,v 



H 



pp ' v + {2( - 1)H? P ' V 



8m 4 p 



'2 

"121-1/2 



2B 



m 



(X - i)w 2 

2m 2 

(2c-d 2 ^ 2B (w-vy +4ol . OI3 



2m 2 p 



+ (2C - 1)H£ 
H ppv +4((l-()H ppv . 



In (|24[1 the following notation has been used: 



B = 



A 2 

T 



16m2C(l-C)' 



(24) 



(25) 



AMPLITUDE OF 77* -> p°p SUBPROCESS 

In this section, we consider the j(q')~f*(q) — > P°(pi)jO°(p2) subprocess. Following [16|, the amplitude of this 
subprocess including the leading twist-2 terms can be written as 

1 



r rll*-^p°p° 



9 =1 n 



g+ u E-(y)V q (y, cos 8, W z ) - ie^E + (y)A q (y, cos 9, W z ) 



where 



E± 



1 i 1 



(26) 



(27) 



1-2/ 2/ 

In (|26[) . the scalar and pseudo-scalar functions (V, A) denote the following contractions 

V q (y, cos0, W 2 ) = e?eiV^H t ^ v (y, C(cos#), W 2 ), 

i 

A q {y, cos9, W 2 ) = J2 eZ4A®Hft A (y, ((cos8),W 2 ). 

i 

The helicity amplitudes are obtained from the usual amplitudes after multiplying by the photon polarization 
vectors 



.4 



- F '(i) F U)rP L ' 



(28) 



Here, in the 77* cm. frame, the photon polarization vectors read 

£ (±) = (0 ±1 zi o^i e (°) = 



191 0,0, 90 



X Q 2 



(29) 



for the real and virtual photons, respectively. 



DIFFERENTIAL CROSS SECTIONS 
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We will now concentrate on the calculation of the differential cross section of . The amplitude of this process 
can be written as 



A 



ee— »eep u p-' 



1,3 



u{l')^u{l) e 



/(<) 



*U) 

s v u{k')Yu{k) 



(30) 



This amplitude depends on the polarization states of the produced p mesons. Due to parity invariance there are 
only three independent sets of helicity (photon) amplitudes which we put to be Ar +t +-\, A( + _^ and A( + Q y Let 
us focus on the leading twist-2 helicity amplitude in the unpolarized electrons case, i.e. the j4( +)+ ) amplitude. 
In this case, the square of the modulus of the amplitude (|30l) can be presented in the "factorized" form : 



|^-ee — >eep°p° | 



u 



\ 2 — \A I 2 

q 



(31) 



and the scattering cross section is 

da. 



d 3 V 



d 3 



Pi 



d 3 



i>2 



d 3 k' 



ee — »eep 1 ' p" 



2S ee (2tt) 3 2Z (27r) 3 2p° (2tt) 3 2 P 9, {2^) 3 2k' Q 



\A e ~{— >ep°p° | ,4 |"4e— >e~f\ 



(32) 



^From now in, we will sum over the polarization states of the p mesons. Separating the differential cross section 
for e7 — > ep°p° subprocess, we are able to rewrite (|32|l in the form : 



d>a ee — >eep° p° 



d 3 V_ X2_, A ,2 , 

/ /4 l-^e^e7l ""Vy— >ep°p° ; 



where 



da„ 



(2n) 3 2l' q 
d 3 p! d 3 p 2 



d 3 k' 



2S ei (27r) 3 2p0 (27r) 3 2p§ (2n) 3 2k' 



■I |-^e7^ep°p°| 



Using the equivalent photon approximation we find the expression for the corresponding cross section 



da, 



ee — >eep u p u 



dQ 2 



dW 2 dcos9 d(j)dx2—Fww{x2) 



da, 



e~f — »ep u p u 



dQ 2 dW 2 dcos0 < 



where the Weizsacker- Williams function Fww is defined as usual as : 

l + {l-x 2 f Q' 2 (x 2 ) l-x 2 
F W w{x 2 ) = z ln- 



2x 2 



3'2 



and the value Q is defined as 



Q' 2 = -a' 2 = 



(1-^2) 



S.-sin a r , 



(33) 



(34) 



(35) 



(36) 



(37) 



The angle ot max in 137(1 is determined by the acceptance of a lepton in the detector (see, for instance, [l2j) and 
the value of the cm. energy of the ee collision y/S ee is 91 GeV at LEP1 and 195 GeV at LEP2. 
In (|35|l . the cross section for the subprocess can be calculated directly ; we have 



da ei _ e pOpQ 
dQ 2 dW 2 dcos6>< 



a 3 1 
16^ S^Q 2 



1 - 



2S ei {Q 2 + W 2 — S ei ) 
(Q 2 + W 2 ) 2 



\A {+ , +) (co S 6,W 2 )\ 



(38) 



where 



A 



(+!+) (cos6»,W^ 2 )| 2 = (|V(cos6»,V^ 2 )| 2 + |A(cos6»,^ 2 )| 2 ) . (39) 



In (|39[) . the squared and p meson polarizations summed functions |V| 2 and |A| 2 read : 
|V(cos0,W 2 )| 2 = ip aiaa (pi)P / ' lft, (p3) 

^V^fd^E^He^iy^cose^^V^J dy 2 E^y 2 )H^ q v (y 2 ,cos8,W 2 ) (40) 

i, 9 J; Q 

|A(cos0,W 2 )| 2 = ip° iaa (pi)P Aft (p2) 



where 



A 

As mentioned before, the expressions (|40|) and l|41|) do not depend on the azimuth <j>. 

The helicity squared amplitude when the integration over cos 9 is implemented may be written as 

l 



F (+,+)( w2 )= J dcos8\A {+ , +) {cos8,W 2 )\ 2 , (43) 
o 

and the cross section takes the form 

i / < M 

x 2 S ee Q 2 7 Q 2 + I^ 2 + Q 2 7 (Q 2 + I^ 2 ) 2 i 1 ' 

In ||33| , the magnitudes of Wl^ and W^ ax are defined by the interval covering most of detected two p° events, 
which in the L3 case at LEP [9j is : 

1.21 GeV 2 < W 2 < 9.0 GeV 2 . (45) 

Since the p meson width is large, the lower limit is a matter of convention, but it can be less than 4to 2 . 
Notice, also, that the integrated function is independent of Q 2 up to logarithms. Besides, the exact Ty 2 

dependence of this quantity remains unknown unless some modeling is used. However, the mean value theorem 
gives the possibility to reduce the three different integrals over W 2 in (|44|l to one integration. Indeed, the mean 
value theorem reads : 

^j#-<?+W / <"*""wn (46) 

W min W min 

W 2 ■ w 2 . 



with two phenomenological parameters (Wx) and (W2). Let us now discuss the status of these parameters. l8 
principle, each of these parameters is a function of Q 2 , but owing to the unknown W 2 dependence of hadronic 
function F/ +>+ \ the Q 2 dependences of given parameters stay out of the exact computations. Therefore we will 
deal with (Wx) and (W2) which will be considered in the sense of average value on the whole interval of Q 2 . 
Notice also that the values of our parameters have actually the same order of magnitude as Q 2 , therefore we 
need to keep (IF)-parameters in the prefactors of and (|4T|) . 

^From H4tjfl and (|47|l . one can see that it is useful to introduce a third phcnomcnogical parameter : 

Cx= J dW 2 l3F {+ , +) {W 2 ). (48) 

W min 

The normalization of our 2pGDA's and consequently the value of our C\ parameter (see, jSQ ) is difficult to 
guess. In the nn case (see, for instance ) , the relation of the second moment of the operator defining the 
GDA to the energy momentum tensor, allowed to relate the value of the second moment of the GDA to the 
total energy carried by quarks in the 7r meson, which is known from the study of parton distributions in the 
meson. This analysis required a modest extrapolation of the meson-pair energy to zero, which was legitimate 
in the two tt meson case with small W, but is very dangerous in the present two p meson case, due to the much 
larger threshold W m i n . Still, one may expect, that some order of magnitude estimate may be provided by such 
an extrapolation, as will be confirmed by the numerical results below. Moreover, such an analysis may give 
indirect access to the parton distributions inside a p meson, and, in particular, to the spin and orbital angular 
momenta carried by quarks. 

Finally, the cross section 144fl is expressed through these three phenomenological parameters in the following 
simple way: 

1 / v 

25a 4 „ f , „ , J 1 



■Cx J dx 2 F ww (x 2 )\ 



dQ 2 36^ " A J vv vv \ x 2 S 2 e Q 2 Q 2 (Q 2 + (W 2 ) 2 ) 2 x 2 S ee Q 2 (Q 2 + {Wx} 2 ) / ' 

v ' 

(49) 

In the L3 Collaboration analysis, the value W belongs to the interval JJS}. Hence we are able to conclude that 
the phenomenological parameters (Wi) and (W 2 ) may take any values inside this interval. 

COMPARISON OF p-MESONS AND LEPTON PAIR PRODUCTIONS 

In this section we compare the p°p° production with the production of a lepton pair /i + p~ in the same kinematics. 
The cross section of the process ej — > ep, + p~ coming through the 77* subprocess is well-known and has the 
form : 

w 2 

TdW 2 (l 1S ei {Q 2 + W 2 -S ei ) \ +u 2 
dQ 2 AS 2 Q 2 J { (Q 2 + W 2 ) 2 ] [ 



r on, (50) 



^e7 

4m?, 



where 



/ ^-(^) = 8(mi±^-^)), ^ = f^- (5i) 



In (|50|) . we can also apply the mean value theorem. But now, by virtue of the known form of the muon functidfi 
1151( 1. the analogues of (Wj) can be explicitly calculated. So, the mean value theorem for the muon case reads : 



I - wrk^ I iW ^"~ (w% 1521 



,00 

4m;, 4m 



Whence, we obtain 



2 

2 



i? 2 (g 2 ), (54) 



(Wi, M ) 2 (Q 2 ) = i?i(g 2 ) - Q 2 , (q 2 + (W 2A , ) ) 2 ^j 
where the following notations are introduced : 

^ = ^W) K = I ^f^(W% L n{ Q 2 ) = J W 2 ^^- (^) 

The Q 2 - dependences of (Wj, (p)) 2 are shown on Fig. 2. The solid line correspond to the function (Wi t ( M )} 2 (<2 2 ) 
and the dashed one to the function {W 2 . (^)) 2 (Q 2 )- The weak Q 2 - dependence of (Wj,^)) 2 justifies the possi- 
bility to use the averaged (W^^) 2 when fitting the data. Besides, a fitting procedure gives us the following 
representation for these functions : 



(Wx <M ) 2 (Q 2 ) = 3.63 - 0.015 Q 2 + 0.471nQ 
(W 2 , M ) 2 (Q 2 ) = 2.90 - 0.024 Q 2 + 0.73 InQ 

Further, the cross section of two p meson production can be written as 



1 wl, . . 



4m? 



where the function M{x2) is defined by the ratio ( one reminds here that 6*67 is proportionnal to the x 2 fraction, 
see ©): 

Zi(as2) 



AA(x 2 ) 



^2 (x 2 )' 



9 ^ Q 2 + (Wi) 2 (Q 2 + (TA7 2 ) 2 ) 2 y /' 

l2fa) = (' ~ ffl+l^W 1 + W 2 + 0^[m)')' ) C ' (5?) 

If one considers the case where Q 2 is large with respect to the invariant mass squared W 2 , we can omit the 
terms of 0(1/Q 2 ) and 0(\nQ 2 / Q A ) in ((57(1 and obtain that M becomes independent of X2 ■ 

25 

TV w — — = 0.01 GeV~ 2 . (58) 
9 /C 



We stress that this value has been obtained providing the value of upper limit W^ ax is fixed (see, Q45fnl 
Actually, the value AT is a function of W^ ax , owing to the W^^-dependence of the integral K. (see, H55[l). and 
this dependence is pretty strong. For example, enhancing W^ ax up to 16.0 GeV 2 the value Af will be halved. 
Thus, the cross section of our process is related to the cross section of lepton pair production as 



d&ee—feep° p° ^ee-tepp 

dQ 5 °- 01 Cl —dQ^- 



(59) 



In the next section we will show that C\ is close to 1.0 GeV . Hence, one can see that the cross section of two 
p meson production is suppressed compared to the cross section of two muon production by a factor which is 
approximately equal to 100. Note that a factor of suppression of the same order was present in the case of two 
7r mesons production jl^ |. 
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Figure 2: The muonic parameters (Wi, ( ll )) 2 (solid) and (W2, (n)) 2 (dashed) as a functions of Q 2 



COMPARISION WITH EXPERIMENTAL DATA 

In the previous section we obtained a simple expression for the two mesons cross section as a function of the 
three parameters (Wi), (W2) and C±. Let us now make a fit of these phenomenological parameters in order to 
get the best description of experimental data. The best values of the parameters can be found by the method of 
least squares, x 2 -method, which flows from the maximum likelihood theorem (see, for instance, As usual, 

the x 2 -sum as a function of parameters is written in the form : 

where P = {(Wi), (W2), Ci} denotes the set of fitted parameters ; a^ xp and af 1 are the experimental mea- 
surements of the cross section and its theoretical estimations ; 5(Ti are the statistical errors. The experimental 
data for the cross section of the exclusive double p° production were taken from the measurement of the L3 
collaboration at LEP Q . 

Minimizing x 2 -sum in (|60[1 with respect to the parameters P we find that the set of solutions P min with their 
confidence intervals are the following : 

C\ = 1.20 ± 0.23 GeV 2 , (Wi) = 3.0 ± 1.9 GeV, (W 2 ) = 1.50 ± 0.09 GeV. (61) 



With this the magnitude of x 2 is equal to 1.40 and, therefore, we have 

A X ff A = - 28 < 1 ( 62 ) 

degree 01 freedom 

The confidence intervals were defined for the case of one-standard deviation. Fig. 3 shows the experimental 
data with the theoretical fit. 
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Figure 3: Cross-section da f , E ^ ef , p o p o/dQ 2 [pb/GeV 2 ] as a function of Q 2 . The theoretical cross-section is plotted for the 
best fitted parameters which are & = 1.2 GeV 2 , (Wi) = 3.0 GeV and (W 2 ) = 1.5 GeV. 

We can see from (|61|) that the confidence interval for the parameter (Wi) covers the whole available interval for 
IF. Moreover the obtained values of (VF2) is quite small. This analysis shows the compatibility of the data with 
the leading twist analysis which we developed. At the same time, one cannot exclude the existence of sizeable 
higher twist contributions to the production amplitude. More theoretical and experimental studies are required 
to make a definite conclusion on higher twist contributions. 



CONCLUSION 

Data thus are in full agreement with the theoretical expectation on the Q 2 behaviour of the cross section. Since 
the effective structure function ^7+,+) is expected to be independent of Q 2 up to logarithms, our analysis of 
the data is a strong indication of the relevance of the partonic description of the process 77* — > p°p° in the 
kinematics of the L3 experiment. 

Much more can be done if detailed experimental results are collected. For instance the angular dependence of 
the final state is a good test of the validity of the asymptotic form of the generalized distribution amplitudes. 
The spin structure of the final state, if elucidated, would allow to disentangle the roles of the nine generalized 
distribution amplitudes. The W 2 behaviour of the cross section may have some interesting features. It depends 
much on the possible resonances which are able to couple to two p mesons. The p + p~ channel may be calculated 
along the same lines. In that case a brehmstrahlung subprocess where the mesons are radiated from the lepton 
line must be added. The charge asymmetry then comes from the interference of the two processes. These items 
will be discussed in a forthcoming publication. More data may be collected at other energies, in particular, in 
e+e" experiments around 10 GeV such as BABAR in SLAC and BELLE in KEK. 
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